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Abstract 

In this paper, we verify Carl de Boor's conjecture on ideal projectors for 
real ideal projectors of type partial derivative by proving that there exists a 
positive ?7 e M such that a real ideal projector of type partial derivative P is 
the pointwise limit of a sequence of Lagrange projectors which are perturbed 
from P up to ?7 in magnitude. Furthermore, we present an algorithm for 
computing the value of such 77 when the range of the Lagrange projectors is 
spanned by the Grobner escalier of their kernels w.r.t. lexicographic order. 
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1. Introduction 

Polynomial interpolation is to construct a polynomial p belonging to a 
finite-dimensional polynomial subspace from a set of data that agrees with a 
given function / at the data set. Univariate polynomial interpolation has a 
well developed theory, while the multivariate one is very problematic since a 
multivariate interpolation polynomial is determined not only by the cardinal 
but also by the geometry of the data set, cf. [H, [ij. 
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As an elegant form of multivariate approximation, ideal interpolation pro- 
vides a natural link between multivariate polynomial interpolation and alge- 
braic geometry jsj. The study of ideal interpolation was initiated by Birkhoff 
[1] and continued by several authors [2I, [s], yl, |6[. 

Actually, ideal interpolation is an ideal projector on polynomial ring 
whose kernel is an ideal. When the kernel of an ideal projector P is the van- 
ishing ideal of certain finite nonempty set S in M'^, P is a Lagrange projector 
on M.[x] := ]R[xi, . . . ,Xd\, the polynomial ring in d variables over M, which 
provides the Lagrange interpolation on H. Obviously, P is finite-dimensional 
since its range is a #S-dimensional subspace of M.[x]. Lagrange projectors 
are standard examples of ideal projectors. 

It is well-known that every univariate ideal projector is an Hermite pro- 
jector, namely it is the pointwise limit of a sequence of Lagrange projectors. 
This inspired Carl de Boorjsf to conjecture that every finite-dimensional lin- 
ear operator on C[a:;] is an ideal projector if and only if it is Hermite. 

However, Boris Shekhtman[3] disproved this conjecture when the dimen- 
sion d > 3. In the same paper, Shekhtman also showed that the conjecture 
is true for bivariate complex projectors with the help of Fogarty Theorem 
(see jsj). Later, using linear algebra tools only, de Boor and Shekhtman 
reproved the same result. Specifically, Shekhtman 10|] completely analyzed 
the bivariate ideal projectors which are onto the space of polynomials of de- 
gree less than n over real or complex field, and verified the conjecture in this 
particular case. 

Let P be an ideal projector that only interpolates a function and its par- 
tial derivatives. Obviously, many classical multivariate interpolation projec- 
tors are examples of P which has applications in many fields of mathematics 



and science, cf.[ll[. Naturally, we wonder whether de Boor' s conjecture is 
true for P or not. 

In this paper, a positive answer is offered to this question by Theorem [2] 
of Section [3] which states that there exists a positive 77 G M such that P is 
the pointwise limit of a sequence of Lagrange projectors which are perturbed 
from P up to ?7 in magnitude, and the proof of the theorem is postponed to 
Sectional the last section of the paper. A further natural question is how to 
determine the value of rj. We propose an algorithm in Section |3]for computing 
the value of such 77 when the range of the Lagrange projectors is spanned by 
the Grobner escalier of their kernels w.r.t. lexicographic order. And then. 
Section 4 is dedicated to some examples to illustrate the algorithm. The next 
section. Section 2, is devoted as a preparation for this paper. 
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2. Preliminaries 



In this section, we will introduce some notation and review some basic 
facts related to ideal projectors. For more details, we refer the reader to 

UMM- 

Throughout the paper, we use No to stand for the monoid of nonnegative 
integers and boldface type for tuples with their entries denoted by the same 
letter with subscripts, for example, o: = (ai, . . . , a^). 

Henceforward, we use < to denote the usual product order on Nq, that is, 
for arbitrary a, /3 G Nq, a < /3 if and only if < z = 1, . . . , d. A finite 
nonempty set () C Nq is called lower if for every q;G0,O</3<q: implies 
/3 G 0. 

A monomial x"" G M\x\ is a power product of the form x^^ ■ ■ ■x^'* with 
cx G Nq. Thus, a polynomial p in M.[x] can be expressed as a linear combina- 
tion of monomials from Supp(p), the support of p, as follows. 



p 



Y,p{cx)x- (1) 



where p{cx) G ]R\{0}. For i G Ng and p G M\x\, if there exists a monomial 
a;" in Supp(p) such that cx' < i, then we denote this fact as p <m *• 

Let P be a finite- dimensional ideal projector on M.[x]. The range and the 
kernel of P are denoted by ranP and kerP respectively. Furthermore, P has 
a dual projector P' on ]R'[a3], the algebraic dual of M.[x], whose range can be 
described as 

ranP' = {A G M'[ic] : kerP C kerA}, 

which is the set of interpolation conditions matched by P. Assume that 
A C W[x] is an M-basis for ranP', then 

kerA := {/ G R[x] : A(/) = 0,V A G A} = kerP. 

We denote by T"^ the monoid of all monomials in M.[x]. For each fixed 
monomial order -< on T*^, a nonzero polynomial / G M.[x] has a unique leading 
monomial LM^(/), which is the -<:-g reatest monomial appearing in / with 



nonzero coefficient. According to [13|, the monomial set 



Ar^(kerA) := {x" G T'^ : LM^(/) f x", V/ G kerA} 

is the Grohner escalier of kerA w.r.t. -<. We denote by ran^P the range of 
P spanned by the Grobner escalier of kerA w.r.t. -<. 
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When P is a Lagrange projector, we have kerA = the vanishing ideal 

of some finite nonempty set S C M'^. In 1995, Cerhenco and Mureddu 14 



proposed an purely combinatorial algorithm named MB for computing the 
Grobner escalier of w.r.t. some lexicographical order on T'^ which is 



denoted by -<iex here. Later, Felszeghy, Rath, and R6nyai|15[] provided a 
faster algorithm, lex game algorithm, by building a rooted tree T(S) of d 
levels from S in the following way: 

• The nodes on each path from the root to a leaf are labeled with the 
coordinates of a point. 

• The root is regarded as the 0-th level with no label, its children are 
labeled with the d-th coordinates of the points, their children with the 
[d — l)-coordinates, and so forth. 

• If two points have same k ending coordinates, then their corresponding 
paths coincide until level k. 

Given finite nonempty point sets E^-^\ C M'^ with = If 

T^E^^^) and T(H*^^)) have same structure, [l^ showed that A/!^j (X(H*^^))) = 
A/'...ms(^))). 

3. Main results 

Let 

6^ : R[x] ^ M : / ^ /(^) 
denote the evaluation functional at the point ^ = ■ ■ ■ ,C,d) £ and let 

Qa Qai^ hOd 



D" : R\x] R\x] : / ^ - — / := — — -/ 



be the differential operator with respect to ck = {ai, . . . ,ad) G Nq with 
D° = I, the identity operator on R[x]. 

Definition 1. Let P be a finite-dimensional ideal projector on R[x]. If 
there exist distinct points . . . , G R^ and their associated lower sets 
D(^),...,0('^) C such that 

ranP' = SpanK{5^(fc) o D° : a G 0^'^^ l<k<fx}, (2) 

namely P only interpolates a function and its partial derivatives, then we 
call P an ideal projector of type partial derivative. 
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As typical examples, Hermite projectors of type total degree and of type 



coordinate degree are both ideal projectors of type partial derivative, cf. [16 



Lemma 1. Let . . . , ^^^^ G M'' be distinct points, and let ()^^\ . . . , O'^'*^ C 
Nq be their associated lower sets. Set 

r]o := min <^ : a G a' G D^^\a ^ a' ,1 < k < I < fi} . (3) 

||q:-q:'||2 J 

Then for arbitrary nonzero h G {—t]o, tjq) C M, t/ie pomt set 

Eh:= {^'^''^ + ha:aed^''\l<k<fi} (4) 

exactly consists of distinct points. 

i=l 

Proof. Suppose that there exist a. G D^'^^ and cx' G D*^'^ with 1 < k < I < n 
such that ^'^'^^ + /la = ^"^'^ + /la' which implies that a ^ a' hy ^^'^^ ^ ^^'^ 
Consequently, we have 

|||W-|(0||, 
\\a. - a'\\2 

which is in direct contradiction to the hypothesis that < \h\ < tiq. □ 

Lemma [1] holds out the possibility of intuitively perturbing an ideal pro- 
jector of type partial derivative to a sequence of Lagrange projectors. 

Definition 2. Let P be an ideal projector of type partial derivative on M.[x] 
with ranP' described by ([2]). For an arbitrary fixed /i G M with < \h\ < r]o 
where tjq is as in ([3]), define to be the Lagrange projector on M.[x] with 

ranP;^ = SpanR{5^(fc)_^.;j„ : a G I < k < fi}. (5) 

Then P^ is called an h-perturbed Lagrange projector of P. 

Remark 1. It is easy to see from ([2]) and ([5]) that 

A := o : a G O^'^^, k = l,...,fi)e (M'[aj])" 

and 

A;, := (<5^(.)+^, : cx G 0^'=), k = l,...,fi)e 
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form M-bases for ranP' and ranP^ respectively, where n = ^^^^ #0^*^-*. 

Moreover, an ordering -<x for the entries of A and A/j will be defined as 
follows: We say 6^(,k) o D" -<a 5^(k') o D"' or 5^(fc)+/,a -<a <^^(fc')+/ia' if 

k < k', or k = k' and o; -< a', 

where -< is an arbitrary monomial order on Nq. 

We are now ready to give one of our main theorem. Theorem |2l which 
states that every ideal projector of type partial derivative on M.[x] is the 
pointwise limit of Lagrange projectors, namely Carl de Boor's conjecture is 
true for this type of ideal projectors. 

Theorem 2. Let P be an ideal projector of type partial derivative on M.[x] 
with lanP' described by and let {Ph, < \h\ < rjo) be a sequence of 
perturbed Lagrange projector of P where rjQ is as in Then the following 
statements hold: 

(i) There exists a positive r/ G M such that 

ranP/j = ranP, VO < < ?7 < ?7o. 

(ii) P is the pointwise limit of Ph,0 < \h\ < 77, as h tends to zero. 

The proof of Theorem |2] will be provided in Section |5l Actually, with 
similar methodology there, we can easily prove the following theorem, which 
is a more general version of Theorem [2J 

Theorem 3. Let P be an ideal projector of type partial derivative from 
C°°(M'^) onto ranP, then there exists Lagrange projector Ph onto ranP such 
that for all f G C°°(]R''), Pf is the limit of Phf as h tends to zero. 

Now, after introducing Definition [3l we have an immediate corollary of 
Theorem |2l 

Definition 3. jsl Let P be an ideal projector from M.[x] onto ranP with 
dimranP = n. Assume that q = {qi, . . . , g„) G Mfa?]" is an M-basis for ranP, 
and the border set dq of q is defined by 

dq:= {1, Xkqu A; = 1, • • • , c?, / = 1, • • • , n} \ {gi, . . . , g„}. 

Then the set of polynomials 

{f-Pf:fe dq} 
forms a border basis for kerP, which is called a q-border basis for kerP. 
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Corollary 4. Let P be an ideal projector of type partial derivative on M.[x], 
and let q be an M.-basis for ranP. Then there exists a Lagrange projector 
onto lanP such that the q-border basis for keiP is the limit of q-border basis 
for keiPh as h tends to zero. 

Theorem [2] tells us that every ideal projector of type partial derivative 
is the pointwise limit of Lagrange projectors. Unfortunately, the converse 
statement is not true in general as the following example illustrates. 

Example 1. Let (P/,, < \h\ < 1) be a sequence of Lagrange projectors 
with 

lanPh = Spanig{l,xi,X2,Xi,XiX2,X2}, 

ranP;^ = Spanig{(5(o,o), ^(o./i), <^(/i,o)5 <^(i,i)5 <^(i,i+/i), (^(i+h,!)}, 

and let P be an ideal projector with 

ranP' = |<5(o,o) o D^^'"), 6^o,o) o 5(o,o) o D^^'^), 

However, {1, xi, ^2, x^, Xia;2, a;2} can not form an M-basis for ranP. Hence, 
(P/j, < \h\ < 1) can not converge pointwise to P, as h tends to zero. 

Consider the bijection 

u : M"^ X N;^ — >{Rx Nof 

(I, a) — > ((^1, ai), . . . , (^d, aa)). 

Let . . . , e R'^ be distinct points and ()^^'^ C be lower sets. 

Then 

n := a) : a G A; = 1, . . . , /i} C (M x No)'' (6) 

is called an algebraic multiset. As mentioned by [13], MB algorithm can be 
applied for the algebraic multiset Q to obtain the Grobner escalier of the 
ideal 

{p e R[x] : 5^(fc) o D"(j9) = 0, a G I < k < fx} 
w.r.t. lexicographic order. 
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Recall Section [21 We have known how to build a d-level tree T(S) from 
a finite nonempty set S G M'^. If the space M"^ is changed to (M x NqY, it is 
easy to see that we can also build a d-level tree T{Q) from algebraic multiset 
Q following the same rules, which makes lex game algorithm involved and 
leads to the following useful lemma. 

Lemma 5. Let P be an ideal projector of type partial derivative with ranP' 
as in ([2]), and let Ph be a perturbed Lagrange projector of P. Let algebraic 
multiset f2 C (M x Nq) be as in ([6D and Eh C R'^ be as in (jl]) . // the rooted 
trees T{Q) and T(S/j) have the same structure, then 

ran^, P = ran^, Ph. 

Next, we can proceed with another main theorem of this paper. 

Theorem 6. Let P be an ideal projector of type partial derivative with ranP' 
as in ([2]), and let {Ph,0 < \h\ < rf) be a sequence of h-perturbed Lagrange 
projectors of P, where rj is obtained through Algorithm^ in the following. If 
the range of Ph is ran^^^^Ph, then the sequence (P/i,0 < \h\ < rj) converges 
pointwise to the ideal projector P , as h tends to zero. 

Algorithm 1. (The range for \h\) 

Input: Distinct points $}'^\ . . . , $}^'^ E R.^ and lower sets . . . , D*^^) C 

K 

Output: A nonnegative number r/ G M or cxo. 

Step 1 Construct algebraic multiset f2 from . . . , ^^^^ and . . . , D*^'^^ 
following ([6]), and then build rooted tree T(f2) from fl in the way introduced 
in Section [21 

Step 2 Suppose that the first level nodes of T{Q) are labeled with the 
points of set £i C M x Nq. 

Step 2.1 If #£i = 1, then ^ oo. 

Step 2.2 If every point in Ci has the same first coordinate or the 
same second coordinate, then 77 ^ oo. 
Step 2.3: Otherwise, set 



rj ■(r- mm 




ii!P,af) and (eJUJ'^) e A 
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step 3 Set ^ 2. 

Step 4 Suppose that the k-th level nodes are labeled respectively with 



the points of sets C 



(1) 



C M X where for each 1 < I < i', the 



nodes labeled with the points in C^^'^ share the same parent. For / = 1, 
and > 2, do the following steps. 

Step 4.1 Set 



rj mm 



Aj) I 



(i) I 



(i) 



7^ a 



(i) 



and {Q_ 



fc+l' '^d-fc+l. 



(i) 



(0 



Step 4.2 If r/' < 77, then rj r/'. 
Step 5 U k = d, then return and stop, 
continue with Step 4. 



Otherwise set k /c + 1, 



Proof. To prove this theorem, by Lemma [5] and Theorem [21 it suffices to 
show that the rooted trees T(S/j),0 < \h\ < rj, and T{fl) have the same 
structure, where is as in and is as in ([S]). Now, with the notation in 
Algorithm [H we will use induction on the number of levels k of the rooted 
tree to prove this. 



When = 1, assume that there exist some 



Oand (eJ\aJ))eA 



such that ^^^^ + haf' 
that/z= \^-^\ 
contradicts 



\a 



a 



The same argument in Lemma [T] shows 



where a^*^ 7^ a^J' and ^''J^' 7^ ^'"j" , which 



d 



\h\ < min 



a 



(i)i 



„(i) 



a«) and 



Hence, the first levels of T(S/j),0 < \h\ < r], and T{Q) have the same struc- 
ture. 

Suppose that the first k — 1 levels of T{Eh),0 < \h\ < r], and T{fl) 
have the same structure. Assume that there exists some 1 < I < and 



iCA^v ^Tu^i) e such that it,^,+haZ,^, 



Ai) 



+ 
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Since (^£^+1, (^i+i, al-fc+i) have common parent, it is 

easy to see that h = l^^l^^^ - ^d-k+i\/\^d-k+i - "i-fc+il ^here a'^^-k+i ^ 



<^d-k 



^■'^ ^ and 7^ ^^Ik+i, which contradicts the fact 

\^d~k+i ^d~k+i\ .Ai) / Aj) U) / U) 

(dVi'^Vi) and 

Therefore, the first k levels of T(S/i),0 < \h\ < r], and 7(17) have the same 
structure. □ 

4. Example 

In this section, we will present several examples to illustrate Theorem El 
Example 2. Assume that is a Lagrange projector with 

ranP;^ = Spanig{5(o,o), \h,o), <^(o,/i), 
Construct the rooted tree of the algebraic multiset 

n = {((0, 0), (0, 0)), ((0, 1), (0, 0)), ((0, 0), (o, i)), 

((1,0),(1,0)),((1,1),(1,0)),((1,0),(1,1))}. 

T{Q) is illustrated in Figure [1] By Algorithm 1, we obtain rj = 1. From 
Theorem IHl we can conclude that {Ph,0 < \h\ < 1) onto Spanjgjl, 2:2, Xi, a;|, 
XiX2,xl} pointwise converges to an Hermite projector P with 

ranP' = {<5(o,o) o D(°'0), ^o) o D^^'"), 5(o,o) o D(°'^), o D(°'°), o D^^'"), 
5(1,1) oD(°'i)}, 

as h tends to zero. 

Example 3. Assume that Ph is a Lagrange projector with 

ranP;^ = Spanig{(5(o,o,o), <^(/i,o,o), <^(o,/i,o), <^(o,o,/i), <^(i,i,i)5 

^{l+h,l,l), S(l,l+h,l), 5{l,l,l+h)}- 
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Figure 1: T{n) of ExampleEl 



Construct the rooted tree of the algebraic multiset 

n = {((0, 0), (0, 0), (0, 0)), ((0, 1), (0, 0), (o, o)), ((o, o), (o, i), (o, o)), 

((0,0), (0,0), (0,1)), ((1,0), (1,0), (1,0)), ((1,1), (1,0), (1,0)), 
((1,0), (1,1), (1,0)), ((1,0), (1,0), (1,1))}. 

T(f2) is illustrated in Figure |2j By Algorithm 1, we compute t] = 1. From 




Figure 2: T(f7) of Example [3] 



TheoremlHl we can conclude that {Ph, < \h\ < 1) onto {1, x^, X2, X2X3, 
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XiX3,x|} pointwise converges to an Hermite projector P with 

ranP' = {5(o,o,o) o D(°'°'°), 5(o,o,o) o D^^'^'^), 5(o,o,o) o D^"'^'"), 5(o,o,o) o D^^'O'^), 
o D(°'°'°), o D(i'°'°), 5(1,1,1) o D(°'i'°), 5(1,1,1) o D(°'°'^)}, 

as h tends to zero. 

Finally, we select test functions 

/l(Xi,X2) = l + (l-Xl)^ + (l-X2)^ 
/2(Xi, Xs, X3) = 1 + (1 - Xif + (1 - + (1 - X3)', 



to illustrate the pointwise convergence of ideal projectors of type partial 
derivative in the above examples. 

For Example 2, when h = 1/10, 1/100, 1/1000, . . ., we have 



100 



1000 



=3 
=3 



385039 
99000 



3439 



■X2 



1000 



Xi 



_ 39984109399 

" 9999000000 
970199 , 



■X2 - 



719 

150" 
3940399 

1000000 



Xl 



-X. 



2i 



249975 
571426287284857 



86 

X\X2 

25 

970199 
165000 



1438 



495 



-X, 



25 



9851 
2500 



XxX2 



3994003999 



-X2 



-x\ + 



997001999 



142857000000000 1000000000 166500000 
998501 142428857 , 



250000 



X1X2 



35714250 



-X 



21 



Pfi =3 - 4x2 - 4xi + 6x2 + 4xiX2 - 4x2. 
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For Example [31 



X . 829 19 19 7 2 o o 80 3 

-T 1 to =4 — : — Xi Xn + 2x9X3 + 2x1X3 H a;q 

w'^ 495 10 10 3 ^ 99 ^ 

„ , , 989299 199 199 37 . 

F 1 /q =4 X3 Xo Xi Xo + 2X9X3 + 2X1X3 

^■'^ 499950 ^ 100 100 33 ^ ^ ^ 

800 3 

^ 9999^^' 

„ , , 998992999 1999 1999 337 , „ 

" 1 /"9 =4 X3 X9 Xi Xo + 2x9X3 + 

499999500 1000 1000 333 ^ 

8000 , 



999999^^' 

-P/2 =4 - 2x3 - 2x2 - 2xi - X3 + 2x2X3 + 2x1X3. 

5. Proof of Theorem [5] 

First of all, we need to relate forward differences of multivariate polyno- 
mials to their partial derivatives. The following formula is quite useful for 
this purpose. 

Lemma 7. Let m G No satisfying i > m > 0. Then 

i:(-i)^(;)(.-.r={o'' :::: m 

Proof. The proof can be completed by induction on m. □ 

Lemma 8. Let ^, /i G M, /i 7^ 0, and i,a & Nq. Then for every monomial x" 
in ]R[x], 

where the remainder 0{h''~^^) is a polynomial in h. 

Proof. From the theory of finite difference(see for example we know 
that 

A^5J{x) = J2i-^y(')s^+hi^-J)f{x) = h% o D7(x) + 0{h'^'), 
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where A is the forward difference operator and /(x) G C*(]R). When f{x) 
is substituted by x"' in this equation, follows immediately. Moreover, 
by Lemma [TJ we can easily check that the remainder 0{h''~^^) in ([8]) is a 
polynomial in h. This completes the proof. □ 

The conclusion of Lemma |8] will be carried over to multivariate 
follows. 

Lemma 9. Suppose that h G M\{0}, | = G R'^, and i = 

(ii, . . . , irf) G Nq. Then for artitrary monomial a;" in M\x\, we have 

2^1 ^) \-\H+Ki^j)^ \ hMi^ oB^x^^ otherwise, ^ ' 

0<j<i ^"^^ S ' 

where (-1)^' = {-lY' ■ ■ ■ {-lY" and (j) = (j^J ■ ■ ■ (j^) provided that j = 
(ji, • • • , id)- 

Proof. First, it follows from Lemma [8] that for every 1 < k < d 



(10) 



Further, we observe that 



\{5^,oW^xl-=5^oWx- (11) 



fc=i 

and 



Finally, we distinguish three cases to prove that the right-hand sides of f ll2p 
and ^ are equal to each other, which will complete the proof. 
Case 1: ex. < i. 

Using f lTOj) and f fTT]) . it is straightforward to verify that 



k=l \jk=0 / k=l 
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Case 2: i ^ CK and a i. 

In this case, there must exist some 1 < k,l < d such that < ik and 
ii < ai. Thus, it is easily checked that 



Case 3: i < ex. 

Let / = max{fc : ik < ak, 1 < ^ < d}. Then, applying (fTOj) and ( ITT]) , we 
deduce that 



ik 



k=l \jk=0 



£,k+h{ik~jk)^k 



I / ik / ■ \ \ / 'fc / • \ 



(,k+h(ik~3k)'-^k 



k=l \jk=0 ^-"^^ / fc=«+l \jfc=0 

I 



^krj.'^k 



fc=l k=l+l 

where the empty product is understood to be 1. □ 

Equation ([9]) makes a connection between the forward difference calculus 
and the differential calculus for multivariate monomials. From Lemma El it 
follows that the remainder 0(/ill*lli"'"^) in ([9]) is a polynomial in h. Equipped 
with these facts, we can establish the relationship between forward differences 
and partial derivatives of multivariate polynomials, which plays an important 
role in the proof of Theorem [21 

Corollary 10. Leti,h,$, be as in Lemma[9l and p G ]R[a;]\{0}. Then 



hMi 2^^ [■y^i+H^-J)P \ S.oD^p, otherwise. ^ ' 

0<j<i ^ ^ 

Proof. Assume that nonzero polynomial p has form (JT]). Since 



i+h{i-3) 



0<j<i a 0<j<i 
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and 

5^ o D'p = ^p{cx)6^ o D^a;", 

ct 

we get 

l^J ^) yjJ''i+Hi-j)P \ h\\'\\^6^oD'p, otherwise, 

which leads to the corollary immediately. □ 

Now, we are ready to prove Theorem [2J 

Proof of Theorem [H We adopt the notation of Definition [1] and Remark 
[TJ Let q = (gi, g2, • • • , Qn) be an M-basis for ranP. Without loss of generality, 
we assume that the entries of A and are ordered ascendingly w.r.t. -<x 
and then denoted as Oi, . . . , o„ and o[, . . . ,o'^ respectively. For convenience, 
we set n X n matrices 

X^q = {oiqj)i<ij<n, A^q = io[qj)i<ij<n, 

and, therefore, for every q G M.[x], n by 1 vectors 

X^q = {oiq)i<i<n, Xlq = (o-g)i<i<„. 

By Corollary [TOl equation (fT3|) can be rewritten as 

I ^ E (-l)^'(i)^S+Mi-i)P' otherwise, 

0<j<i 

which implies that for fixed 1 <k < fi and i G 0'^'^\ 5^{k) o D*p can be linearly 
expressed by {5^(fe)+/jjP : I G 0*^*^^} U {0{h)} since D^^^ is lower, and moreover, 
the linear combination coefficient of each S^{k)_^_fiiP is independent of p G M.[x]. 
Thus, it turns out that there exists a nonsingular matrix Tp of order n such 
that 

A^jAfg] := T, [A^glA^g] = [A^q|A^g] + [E,\e,] , (14) 

where each entry of [-E/i|e/i] is either oi 0{h). As a consequence, the linear 
systems 

{Kq) X = Kl and [Xlq] x = X^q 
16 



are equivalent, namely they have the same set of solutions. 

(i) From (fT4|) . it follows that each entry of matrix X[q converges to its 
corresponding entry of matrix X^q as h tends to zero, which implies that 

linidet (^^q) = det (X^q) . 

Since det(A^q) 7^ 0, there exists 77 > such that 

det (^A^) ^0, 0<\h\< rj. 

Notice that f[T^ directly leads to rank (^Xj^q^ = rank (A^q) , 

lanPh = Spauj^qr, < \h\ < r], 

follows, i.e., q forms an M-basis for ranP/j. Since q is also a basis for ranP, 
we have 

ranP = ranP/j, < \h\ < rj. 

(ii) Suppose that Xh. and x be the unique solutions of nonsingular linear 
systems 

{Xlq)x = Xlq (15) 

and 

{X^q)x = X^q (16) 
respectively, where < < 77. It is easy to see that 

Phq = qxh and Pq = qx. 

Remark that, as /i — )■ 0, P is the pointwise limit of Ph if and only if Pq is 
the coefficientwise limit of Phq for all q G M.[x]. Therefore, it is sufficient to 
show that for every q G M.[x], the solution vector of system (fTS!) converges to 
the one of system fll6p when h tends to zero, namely 

lim Xh = X. 

By (fT4l) . the linear system 

■Xlq)x = Xlq (17) 



17 



can be rewritten as 



[X^q + Eh)x= {X^q + Eh) . 



Since system (fT7|) is equivalent to system (fT5|) . Xh is also the unique solution 
of it. Consequently, applying the perturbation analysis of the sensitivity of 
linear systems (see for example [18|, p.SOff), we have 



\xh - x\\ < 



\eh- Ehx\\+0{h^ 



Since each entry of vector eh 
lim II — x\\ 

the theorem. 



Eh.x is either or 0{h), it follows that 
0, or, equivalently, lim x^ = x, which completes the proof of 

/i-*-0 

□ 
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